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ABSTRACT 


The  spectrum  of  seismic  body  phases  radiated  from  sources 
within  tho  earth's  crust  is  distorted  by  multiple  re¬ 
flections  in  the  crustal  layers.  To  estimate  the  crustal 
distortion  crustal  transfer  functions  for  P-v/aves  have 
been  computed  for  sources  having  a  simple  radiation 
pattern  in  the  infinite  medium!  the  buried  explosive 
point  source  standing  for  an  underground  explosion,  and 
the  vertical  point  source  acting  at  the  free  surface  of 
the  crust  and  standing  for  an  atmospheric  explosion.  The 
crustal  model  is  taken  from  Finland. 

Radiation  patterns  have  been  computed  for  P-waves  in  the 
frequency  range  from  0  to  0.5  cps.  An  attempt  ha3  been 
made  to  separate  the  effect  of  the  free  surface  from  that 
of  *-he  layered  crust  in  case  of  the  explosive  source. 

The  latter  effect  starts  to  be  noticeable  at  frequencies 
above  0.3  cps.  Using  tne  model  for  the  upper  mantle  model 
in  Europe  (Lehmann,  1959)  the  radiation  patterns  have  been 
ins  formed  into  amplitude-distance  functions.  Implications 
.  the  interpretation  of  experimental  amplitude-distance 
curves  of  P-waves  are  discussed. 
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1.  INTRODUCTION 


For  the  study  of  the  nonelastic  behavior  of  the  upper 
mantle  we  want  to  use  body  phases  from  explosive  sources 
since  their  radiation  pattern  is  comparatively  simple. 
These  explosive  sources  are  located  within  or  on  top  of 
tne  crust,  Multiple  reflections  within  the  layered  earth's 
crust  are  distorting  the  primary  signal,  originally 
radiated  by  the  source.  Therefore,  the  signal  entering 
the  upper  mantle  r.t  the  base  of  the  "source"  crust  is 
different  from  the  primary  signal.  Its  spoctru u  will 
depend  on  the  elastic  parameters  and  the  density  of  the 
layers  within  the  crust,  its  thickness,  the  anale  of 
incidence  and  the  depth  of  the  source.  Also  most  earth¬ 
quakes  occur  within  the  crust  of  the  earth.  As  an  example 
the  depth  distribution  of  hypocenters  in  South-Western 
Germany  is  given  in  Figure  1  (see  Schneider,  1964)  , 
Therefore,  we  have  to  ask  to  what  extent  will  the  ra¬ 
diation  pattern  of  shallow  seismic  sources  be  effected 
by  the  source  crust.  If  necessary  corrections  for  the 
effect  of  the  source  crust  have  to  be  applied. 

v. 

The  vertical  point  source  acting  on  the  free  surface 
of  the  layered  earth's  crust  will  be  used  to  estimate 
the  influence  of  the  crust  on  the  spectrum  of  seismic 
body  waves  radiated  by  atmospheric  explosions.  The 
explosive  point  source  buried  at  some  depth  within  the 
crust  is  our  model  for  underground  explosions.  After 
computing  the  radiation  pattern  for  the  far-field  at 
different  frequencies  we  will  investigate  how  the  normal 
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Depth  distribution  of  earthquakes  in  S.h.  Gerrany:  nur.iber  of  hypocenters  per  kilorr.eter 
of  depth.  (Data  fron  Schneider,  1964) 


amplitude-distance  curve  is  effected  by  the  distorted 
radiation  pattern. 

The  formulas,  for  the  transfer  function  for  P-waves  for 
a  system  consisting  of  a  point  source  within  a  layered 
medium  over  a  homogeneous  half-space  have  been  derived 
by  Fuchs  (1966)  using  the  theory  of  Ilarkrider  (1964) 
for  wave  propagation  from  point  sources  in  layered 
media.  Those  formulas  arc  valid  for  the  far-field  within 
the  half-space.  The  following  assumptions  regarding  the 
model  "point  source  in  the  crust"  have  been  madet 

(1)  The  actual  source  is  approximated  by  a  point  source; 
i.e.  the  dimensions  of  the  source  region  con  be 
taken  as  small  compared  to  the  distance  to  the 
next  boundary.  With  this  limitation  in  mind  the 
point  source  may  be  placed  at  any  depth  in  a 
layered  medium, 

(2)  The  cru3t  of  the  earth  is  approximated  by  a  system 
of  (n-1)  homogeneous  isotropic  layers  with  plane 
parallel  interfaces.  The  top  layer  is  specified 

to  be  a  free  surface.  The  angle  f  is  the  angle  of 
incidence  of  P -waves  in  the  mantle. 

*  th 

(3)  The  mantle f  the  n  layer,  will  be  regarded  to  be 
"relatively  homogeneous".  By  this  we  mean,  generally 
speaking,  that  the  elastic  properties  of  the  mantle 
vary  more  gradually  than  they  do  in  the  crust. 
Therefore,  we  choose  the  n^  layer  to  be  a  homo¬ 
geneous  half-space. 

(4)  The  solution  becomes  invalid  for  f  ■  90°  duo  to 
interference  with  the  interface  wave  guided  at  the 
base  of  the  crust. 


The  vertical  noint  force  at  the  free  surface  of  the 
layered  crust  is  taken  as  a  first  approximation  of  an 
atmospnoric  explosion.  This  point  force  is  generating 
a  Fourier  transformed  dilatational  displacement  potential 
1  (r,z)  at  a  point  (r,z)  in  the  half-space* 


❖ (r, s) 


cos  y  • 


(1) 


This  formula  is  taken  from  Fuchs  (1966,  equation  (2,81)) 
'./it h  slight  modifications. 


The  following  notation  is  used: 

L(u)  ■  Fourier  transformed  source  function 

k  ®  "  dilatational  wave  number  in  the 

an  a  ,  ,  _ 

n  half-space 

u  ■  2nv  «  angular  frequency 

Y  «=  angle  of  incidence  in  the  half-space 

R  a  distance  from  the  plumb  point  of 

the  source  at  the  base  of  toe  crust 
to  the  receiver  at  (r,z). 


The  quantity  hv  contains  the  effect  of  the  layered  crust. 

It  is  a  function  of  frequency,  elastic  parameters,  thickness 
and  density  of  the  layers  and  the  angle  of  incidence  y 
into  the  half-space.  It  is  computed  from  the  crustal 
matrix  as  given  by  Fuchs  (1966)  , 
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From  equation  (1)  the  I  uricr  transformed  total  uila- 
tational  displacement  u^  (r,z)  at  (r,z)  is  deriveui 


u ( r ,  z )  =»  .  cos  y  •  a 

1  71  \ 


-ik  .R 

an 


(2) 


The  factors  T  ( wt  y)  =  ■*—  cosy  »  A,,  (3) 

may  be  regarded  as  the  transfer  function  for  the  system 
vortical  point  source  at  the  free  surface  of  the  layc?rd 
crust.  T(w/y)  has  been  computed  on  the  13: 1  7094  of  the 
German  Computer  Center  at  Darmstadt .  We  have  chosen  the 
model  iJU?v  1  from  Finland  listed  in  Table  1  (Penttilu, 

1965)  as  an  estimate  for  the  crust  under  Wovaja  Scmlja. 


Depth 

(km) 

Thickness 

(jam) 

(Lm/s) 

3 

(kn/s) 

Density 
(g/err,  ) 

2.2 

5,73 

3.52 

2 , 30 

20.3 

29.0 

13.6 

5.95 

3.52 

2,3ft 

3.2 

6.37 
o  1 23 

3.72 

4.50 

2.95 

3.43 

Ta^lo  1 : 

Crustal  ho del 

UUR  1 

-  Finland 

(Penttilu,  1965) 

riguro  2  depicts  T(v,y)  for  various  angle  of  incidence 
(y  53  20°,  30°,  40°  ar.u  50°)  in  the  frequency  range 
O^v—  0,5  cps.  Strong  deformation  of  the  radiation 
pattern  is  to  be  expected  above  0.1  cps.  The  maxima  of 
T(v  ,  20°)  are  nearly  coinciding  with  the  minima  of 
T(v  ,  50°)  and  vice  versa. 

The  radiation  pattern  for  tne  frequencies  V  =  0  to  0.5  cps 
are  given  in  Figure  3  at  increments  of  0,05  cps.  Only 
the  right  half  sides  of  tne  radiation  patterns  have  been 


ction  T(v,y)  for  t'.ac  vertical  point  source  on  top  of  crust 
°,  40°  and  5 0°) 
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of  crust  hU-l  1  (see  Table  1)  (v  «*  0,0  to  0,5  cps)  , 
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drawn,  The  left  half  sides  arc  symmetrical,  At  0.0  cps 
the  radiation  pattern  corresponds  to  that  of  tho  vortical 
point  force  acting  at  tho  free  surface  of  a  homogeneous 
half-space  of  subcrustal  material.  As  tho  frequency  is 
increasing  tho  radiation  pattern  is  gradually  distorted. 

At  some  frequencies  (o.g.  0.10,  0.2 0,  0.30  cps)  tho 
Maximal  energy  is  no  lonjer  radiated  at  y  =  0°,  that  is 
vertically  into  the  half-space,  hut  at  angles  y  between 
20°  and  *10°.  The  effect  on  tnc  amplitude-distance  function 
will  he  discussed  in  section  4. 


3.  Tab  TAAkSFLA  FollCTIOH  FO 11  Ah’  EXPLOSIVE  POIh’T  SOURCE 


AT  .TPTii  IF  THE  CRUST 


The  transfer  function  for  an  explosive  point  source 
buried  at  some  depth  in  the  crust  has  also  been  derived 
by  Fuchs  (1966).  This  transfer  function  T^c^y)  is  defined 
as  tne  ratio  of  the  displacement  caused  by  the  explosive 
point  source  in  the  crust  to  the  displacement  that  would 
ha-'e  been  caused  if  the  same  source  had  been  located  in 
an  infinite  homogeneous  medium  of  subcrustal  material  at 
tne  plumb  point  of  the  source  at  the  base  of  the  crust. 

It  is  again  a  function  of  frequency,  elastic  parameters, 
thickness  and  density  of  the  crustal  layers  and  the  angle 
of  incidence  in  the  half-space. 

As  an  example  of  an  underground  explosion  the  transfer- 
function  for  an  explosive  point  source  at  a  depth  of 
500  u  within  the  crustal  model  EUR  1  has  been  computed. 

Tho  rauiation  pattern  for  the  frequencies  v  =*  0,0  to 
0.5  cps  at  0,05  cps  increments  are  given  in  Figure  4. 


hndiation  patterns  for  the  buried  explosive  point  source 

in  one  crust  h"Jh  1  (Table  1)  at  depth  jOO  n  (v  «=  J,0  to  u.3  cps) 
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The  radiation  pattern  of  the  explosive  source  differs 
significantly  from  that  of  the  vertical  point  force. 

In  the  case  of  the  buried  explosive  source  a  strong 
primary  reflection  is  generated  at  the  free  surface 
of  tne  crust.  Its  interference  with  the  primary  wave 
is  dominating  the  radiation  pattern.  This  strong  re¬ 
flection  is  absent  in  the  case  of  the  vertical  point 
source  at  the  free  surface. 

At  0,0  cps,  where  the  crust  is  ignored,  the  radiation 
pattern  i  s  that  of  an  explosive  dipo.le  in  an  infinite 
medium,  voice  radiates  no  energy  in  the  direction  of 
tne  dinole  axis.  As  the  frequency  increases  the  finite 
eistance  of  the  source  from  the  free  surface  allows 
energy  also  to  be  emitted  at  y  =  0°,  The  direction  of 
the  radiation  of  maximal  energy  is  a  function  of  frequency • 

It  is  important  to  got  an  estimate  of  that  part  of  the 
transfer  function  which  is  caused  by  the  free  surface 
alone  and  that  part  which  is  caused  by  the  layering  of 
tne  crust.  While  the  first  part  is  the  s^raa  for  all 
crustal  models  -  provided  the  source  is  located  at  a 
corresponding  depth  -  the  second  part  will  vary  with 
the  cr  .stal  model,  Tor  comparison  we  have  placed  the 
explosive  point  source  in  a  homogeneous  half-space  of 
tne  s uterus tal  material  of  tne  UuA  1  model,  but  at  a 
dcpt.i  of  715  m.  This  depth  has  been  chosen  in  such  a 
way  that  the  travel  time  for  the  P-v/ave  from  the  source 
to  the  free  surface  remains  the  same  at  normal  incidence 
as  in  the  case  of  the  source  in  the  h’b'P  1  crust  at  a 
depth  of  500  m.  The  radiation  patterns  are  drawn  in 
Figure  5.  Especially  for  frequencies  larger  than  0,3  cps 
the  radiation  patterns  for  the  explosive  source  buried 
in  the  homogeneous  half-space  are  much  smoother  than 
for  the  same  source  buried  in  the  layered  crust. 


;  x-r.  j:  Radiation  patterns  for  t.\-~  juried  ^xolosiva  point  source  in 

t.ic  iioir.ofjcneouc  half-space  at  ueptl  15  m  (v  a  0.0  to  j,5  ops) 
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i'o  study  t.io  effect  of  t.ic  source  depth  the  explosive 
source  has  also  .jean  placed  at  a  depth  of  25  kin  in 
tne  crust  h’Uh  1.  The  redaation  patterns  are  given  in 
Figure  G.  how  strong  energy  is  radiated  vertically 
into  tne  naif-space  (y  =  0°)  at  comparatively  low 
frequencies  of  about  0,05  cps,  here  the  direct  wave 
from  the  source  and  the  reflection  from  the  free 
surface  arc  interfering  constructively,  To  estimate 
the  effect  of  the  free  surface  alone#  the  radiation 
pattern  of  the  explosive  source  in  a  homogeneous 
half-space  of  subcrustal  material  is  depicted  in  Figure  7, 
Tiie  depth  of  the  source  is  32,4  km  with  equal  travel  time 
to  the  surface  at  normal  incidence, 

h  comparison  of  the  radiation  patterns  of  Figure  G 
and  7  shows  tnat  there  arc  considerable  differences. 

Only  in  the  range  from  0  to  0,1  cps  the  patterns  show 
similar  shape.  For  larger  frequencies  the  influence 
of  tiie  crust  causes  appreciable  deformation  of  the 
radiation  patterns. 


4,  Tdd  1. 1  FLU;  dICd  OF  T1IIJ  PADIATIOU  PATTFili  I  OK  Tub 


A’  i?  L  If  UD  I>  L>  I S TAdCE  FUNCTION 


Aadi  »tion  patterns  cannot  be  observed  directly,  What 
we  do  observe  arc  amplitudes  as  a  function  of  cpiccntral 
ui stances,  fo  understand  how  tne  amplitude  distance 
function  is  effected  by  the  radiation  pattern#  especially 
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for  waves  penetrating  the  upper  mantle,  v/e  have  con¬ 
verted  the  angle  of  incidence  y  at  the  base  of  the 
crust  into  cpiccntral  distances  A.  The  best  model  known 
for  the  upper  mantle  in  Europe  is  the  one  derived  by 
Lehmann  (1959) ,  It  is  listed  in  Table  2  and  depicted 
in  Figure  3.  The  same  figure  also  gives  the  uranchcs 
of  the  T  -  A  -  diagram  of  the  refracted  and  overcritically 
reflected  P-waves  for  this  model*  Travel  times  have  been 
reduced  with  a/6 ,0(krn/ (km/sec) )  .  The  branch  A  corresponds 
to  the  waves  refracted  between  the  li-discontinuity  at 
a  depth  of  35  km  and  the  first  order  discontinuity  at 
220  km.  branch  b  is  the  ovcrcritical  reflection  from 
this  discontinuity  at  a  depth  of  220  km,  branch  C 
corresponds  to  rays  penetrating  tne  upper  mantle  at 
depths  larger  than  220  kns. 

The  y  -  A  -  curve  for  this  model  of  the  European  upper 
mantle  has  been  computed  for  a  surface  focus  and  is 
shown  in  Figure  9.  The  oranchos  have  again  been  labelled 
A, b  and  C  using  the  notation  of  Figure  7, 

This  y  -  A  -  curve  has  now  been  used  to  convert  the 
radiation  patterns  T(y,v)  into  an  amplitude-distance 
function  T(A,v)  of  cpiccntral  distance  A  and  frequency  v. 

T  ( A , v )  docs  not  contain  a  correction  for  geometrical 
s  v-'^ading  and  amplitude  splitting  at  first  order  dis¬ 
continuities.  It  is  so  to  speak  a  projection  of  the 
radiation  pattern  on  the  free  surface  of  tne  spherical 
carta  following  the  ray  paths  in  the  upper  mantle. 

Tiie  amplitude-distance  function  T(A,v)  has  been  plottec. 
in  Figure  10  for  tne  explosive  point  source  at  500  m 
depth  in  the  crustal  model  EUR  1  and  for  the  vertical 
point  force  acting  on  the  sane  crust  in  Figure  11,  Tne 
frequency  range  is  0,0  to  0,5  cps  in  increments  of  0,05  cps. 
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(km)  (km)  (km/sec) 


0 

G370 

6.2 

35 

6  335 

6.2 

35 

G335 

9.0 

55 

6315 

3.12 

95 

6275 

3.12 

159 

6211 

3.12 

220 

6150 

3.12 

220 

6150 

3.40 

236 

6034 

3.63 

3-19 

6021 

3 . 35 

412 

5953 

9.24 

476 

5394 

9.54 

539 

5331 

9.36 

G02 

576-3 

10.13 
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Tnole  2: 


P-./ava  velocities  of  the  upper  mantle  in  Europe 
(Lehmann,  1959) 
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ri j.Q:  ."urjie  of  incidence  y  at  base  of  crust  to  distance  A  relation 
for  European  upper  mantle  model  (see  i igurc  3) « 


Fig .  10:  hnplituue-distance  function  T(A,v)  fro a  radiation  patterns  of  explosive  point  source  at  depth  500  :r. 
in  crust  :.TUR  1  (compare  Figure  4) 


2 


Amplitude -distance  function  T(A#v)  from  radiation  patterns  of  vertical  point  source  on  top  of 
crust  AUA  1  (compare  Figure  3)  « 
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Tug  following  points  deserve  special  attention  in  the 
case  of  the  explosive  source  (Figure  9) t 

1. )  Within  the  frequency  range  0  to  0.5  cps  the 

radiation  pattern  of  the  explosive  source  is 
causing  a  relative  increase  of  amplitudes  with 
distance  along  the  refracted  branch  A,  This 
implies  that  any  observed  amplitude  attenuation 
along  branch  A  should  be  even  stronger  if  corrected 
for  the  effect  of  the  radiation  pattern. 

2. )  Amplitudes  along  the  .branch  3  are  slightly  de¬ 

creased  with  distance  by  the  radiation  pattern, 

3. )  The  most  significant  effect  of  the  radiation  pattern 

is  to  be  observed  along  the  branch  C.  The  maximum 
of  the  radiation  pattern  occurs  on  this  branch 
at  cnicentral  distances  A  between  11°  and  21° , 

Fcr  most  frequencies  amplitudes  are  strongly 
decreasing  beyond  this  maximum.  This  is  especially 
true  at  tiie  following  frequencies;  0,0,  0,05,  0.1, 
0,2,  0,25,  0,35,  0,45  and  0,5  cps.  Only  at  0,4  cps 
amplitudes  arc  increasing  with  distance.  In  the 
case  of  0,0,  0,05,  0,2  and  0,25  cps  the  amplitude 
decrease  amounts  to  about  40  »  between  distances 
A  =  11°  and  27°, 

It  should  bo  noted  that  the  distance  range  of  this 

strong  decrease  of  amplitudes  coincides  with  a  minimum 

in  experimental  amplitude  distance  curves  for  P  and 

n 

P-wavcs  oy  Gutenberg,  Aichter  (1956),  de  Bremaeckcr 
(1955)  and  Vane);,  Stclzner  (1962)  ,  It  is  possible  that 
this  minimum  is  not  caused  oy  geometrical  spreading  of 
t.ic  ray  bundle  or  non-elastic  behavior  of  the  penetrated 
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region  but  simply  by  tae  radiation  pattern  of  the 
source  in  the  crust.  The  amplitude-distance  function 
for  r -wave s  radiated  from  the  vertical  point  source 
(Figure  10)  differs  considerably  from  that  of  the 
explosive  source-  While  the  behavior  of  the  amplitudes 
along  branches  A, 3  has  the  same  tendency  as  in  the 
case  of  tiie  explosive  point  source  -  an  apparent 
increase  of  amplitudes  with  distance  caused  by  the 
radiation  pattern  -  the  amplitudes  along  branch  C 
display  a  strong  increase  of  about  100  t  in  the 
distance  range  L  from  11°  to  27°,  Therefore,  experiment¬ 
ally  determined  amplitude-distance  function  from  under¬ 
ground  explosions  should  show  stronger  attenuation  than 
corresponding  functions  from  atmospheric  explosions. 

We  are  planning  to  compare  experimental  data  from  both 
typ  .  j  of  explosions  to  study  the  theoretical  predicted 
effect,  'we  shall  also  investigate  how  variations  of  the 
crustal  parameters  will  influence  the  amplitude  distance 
function. 
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”  D  I 


Administrative  Inrr.rnation 


a)  rachnical  status 


equipment  for  tne  digitization  of  paper  and  film 
seismograms  has  been  installed  at  our  institute  during 
December  ’966,  At  pre  ent  we  arc  testing  the  data  on 
the  punched  paper  tape  and  the  transfer  to  our  programs 
for  spectral  analysis  at  the  Ib.'i  7094  of  the  German 
Computer  Cento.'  in  Darmstadt.  The  general  overhaul 
of  tiie  various  short-period  seismometers  at  the 
Graofonborg  Observatory  (GGGR)  has  continued. 

«)  hajor  Accomplishments 

Our  efforts  concentrated  on  a  study  of  the  distortion 
of  the  spectrum  of  seismic  body  waves  by  the  crust 
in  the  neighborhood  of  the  source.  Radiation  patterns 
and  amplitude -distance  functions  have  been  computed 
for  the  buried  explosive  source  and  the  vertical  point 
source  at  the  free  surface  of  the  crust. 

The  routine  maintenance  at  the  technical  installations 
of  the  Graefenberg  Observatory  continued.  The  systematic 
measurements  of  seismic  background  noise  at  various 
places  in  Germany  are  still  in  progress, 

c)  Problems  encountered 


fhe  GbOTbCii  low-speed  tape  recorders  are  still  out  of 
order.  They  will  be  returned  to  Dallas  (Texas)  for  repair. 
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d)  Actions  required  by  AFQSA  and  E OAR 

■'o  immediate  action  of  AFOSR  and  GOAR  is  required 
at  this  time, 

e)  Future  olans 

work  v;ill  continue  along  the  lines  of  our  original 
project  outline.  The  special  study  cf  the  offoct  of 
radiation  patterns  on  amplitude  distance  functions 
will  be  carried  on, 

f )  Personal  and  Organization 

Since  September  1#  1966  hrs,  Uieisch  is  employed  as 
a  half-day  secrotary  at  the  Gracfenberg  Oesorvatory, 
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»1  ABSTRACT 
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..to  spectrum  of  seismic  body  phases  radiated  from  sources  within  tite 
carta's  crust  is  distorted  by  multiple  reflections  in  the  crustal  layc.s 
Vo  estimate  the  crustal  distortion  crustal  transfer  functions  for 
i’-waves  have  ween  computed  for  sources  having  a  simple  radiation  pattern 
in  tno  infinite  medium:  the  ouried  explosive  point  source  stanaing  for 
an  underground  explosion,  and  the  vertical  point  source  acting  at  tne 
free  surface  of  the  crust  and  standing  for  an  atmospheric  explosion, 

Fho  crustal  model  is  taken  from  Finland, 

Radiation  patterns  have  been  computed  for  P-waves  in  the  frequency  range 
from  0  to  0,5  cps ,  An  attempt  iias  been  made  to  separate  the  effect  of 
tne  free  surface  from  that  of  the  layered  crust  in  case  of  the  explosive 
scarce.  The  latter  effect  starts  to  oc  noticeable  at  frequencies  aeove 
0,5  cps.  Using  the  model  for  tne  upper  mantle  model  in  Europe  (Lehmann, 
1959)  tne  radiation  patterns  have  been  transformed  into  amplitude- 
eistanco  functions.  Implications  for  the  interpretation  of  experimental 
amnlitudo-distancc  curves  of  P-waves  arc  discussed. 
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